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Abstract. 

A recent reanalysis of world proton and deuteron structure function measurements showed that 
^vq a significant amount of the apparent model dependence in the extraction of the neutron structure 

function was related to inconsistencies between the kinematics of the data and those assumed for 
the calculation, suggesting that the true model dependence is smaller than commonly believed. We 
r-r present a detailed comparison of the neutron structure function as extracted using different models, 

with care taken to ensure that all other aspects of the comparison are done consistently. The neutron 
structure function is extracted using a fit to these data evaluated at fixed Qq=16 GeV 2 . We compare 
the results obtained using a variety of N-N potentials and deuteron binding models to determine the 
model dependence of the extraction. As in the recent extraction, F^nlF-ip falls with x with no sign of 
D plateau and follows the low edge of the wide range of earlier Fm extractions. The model-dependent 

1 { uncertainty in F2 n /F2 P is shown to be considerably smaller than previously believed, particularly at 

Q large-x. 
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v£ INTRODUCTION 

O 

Because of the experimental impracticality of performing precision scattering experi- 
ments with a neutron target, the determination of the structure functions of the neutron 
requires 1) inclusive DIS scattering data from both proton and deuteron targets and 2) a 
theoretical model that convincingly describes the impact on the structure function of the 
proton and neutron when bound together. Previous extractions of the neutron structure 
function, and more relevantly the ratio of the structure function of the neutron to that 
of the proton, F^nl^ip, were performed using different subsets of the world data on Fi p 
and F2d with varying treatments of binding and kinematics. These produced significant 
variation in the resulting ratio, particularly at large-*. Some of these are shown in Figure 
1 . Because these different approaches yielded very different results, it was concluded 
that the model dependence of such an analysis was so large that no reliable information 
on the neutron could be obtained at large -x. However, as discussed in [1], a significant 
amount of the variation between results is related to self-inconsistencies in the extrac- 
tions, giving the impression that our knowledge of this quantity is significantly worse 
than it will actually be shown to be. By carefully and uniformly treating the source data 
and restricting deuteron models to a set that is modern and compelling, a reasonable and 
significantly smaller uncertainty estimate can be made. 
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FIGURE 1. Several prior extractions of the ratio of the neutron to proton structure functions showing 
a large range, particularly at high-x. The models shown are: off-shell (black squares) and on-shell (solid 
red diamonds) extractions from Ref. [2], pure Fermi-motion (open red diamonds) [3], and nuclear density 
dependent EMC effect based models (solid and open blue circles) from Ref. [3] and [4] respectively. The 
arrows on the right-hand side of the figure indicate various theoretical predictions of this ratio in the x = 1 
limit [2]. Exact spin-flavor SU(6) symmetry implies a value of 2/3, symmetry breaking through vector 
(5=1) diquark suppression at large-x gives a value of 1/4, and a value of 3/7 is obtained if only the 
z-component, S z = 1 of the vector diquark's spin determines its suppression. 



EXTRACTING F 2 n/F 2p 

We begin with a summary of the extraction performed in Ref [1], as this provides the 
starting point for this comparison. The proton input, F2 P (x, Q ), used in this analysis 
comes from a parameterization provided by M. E. Christy, which is fit to a large body 
of experimental experimental results. The deuteron to proton structure function ratio, 
Fid/Fjp, is taken from SLAC, BCDMS, and NMC measurements covering a range in 
Q 2 from 3 to 230 GeV 2 . Data points for which 6 < Q 2 < 30 GeV 2 and W 2 >3.5 GeV 2 
are interpolated to a common of 12 GeV 2 (typically a 0.5% correction). The data 
points and fit function are shown in the left panel of Figure 2. 

The nuclear modification of the proton and neutron was computed in a light cone 
impulse approximation, without any assumption of scaling or neglecting k±, using the 
CD-Bonn [5] effective two-body potential. This model is applied to the proton data and 
a R np (= F2n/F2 p ) ratio of the form 

R np (^Q 2 ) = ( Pl +p 2 $)+P3 e-P* « + P5 e~" ^+p 7 max(0,^ -p 8 ) 2 , 

where % = 2x/(l + yi + 4M 2 x 2 / Q 2 ) is chosen as the scaling variable to reduce the Q 2 
dependence. The free parameters are tuned using MINUI T to best reproduce the F^ data 
in conjunction with the parameterized proton structure function. The result of this proce- 
dure is the curve in the right panel of Figure 2. Once R np is determined, smearing ratios 
S p = F2 p /F2p and S n = Fin/Fin* where the tilde indicates the structure function of the 
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FIGURE 2. Left panel: F^/Fip data from the range 6 < Q 2 < 30 GeV interpolated to a common 
2q=12 GeV 2 shown with a parameterization. Right panel: The extracted parameterization of FinjF%p 
with normalization uncertainty bands propagated from the proton and deuteron data sets. The data 
points are generated from the original Fin/Fip values and the smearing ratios produced using the R np 
parameterization. Figures reproduced from [1]. 



bound nucleon, are trivial to determine. Using these ratios, the individual Rd p — F%^jF2p 
points can be converted to R np values, along with their individual uncertainties. In ad- 
dition, detailed estimates of the systematic uncertainties associated with the extraction 
procedure and data normalization have been parameterized. The individual R np points 
as well as the systematic uncertainties are shown as the points and dotted bands in the 
right panel of Figure 2. 
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FIGURE 3. The resulting parameterization (red line) and data points (solid red circles) compared 
with previous calculations. An off-shell model [2] (open green squares), on-shell model [6] (solid black 
squares), and a modified EMC effect model (open blue circles) [3] are shown. Reproduced from [1]. 



Figure 3 compares the results of this extraction to previous results, with the new data 
points rebinned to approximately match the previous binning. These results fall on the 




FIGURE 4. Left panel: The calculated ratio Fq^ jFip evaluated at the average Q 2 for the given x value 
(red), and at a common Qq = 12 GeV 2 (black). Right panel: The proton smearing ratio, S p for values of 
Q 2 ranging from 4.7 to 23.6 GeV 2 (corresponding to the average Q 2 value of each of the Rj p points). The 
dark central line corresponds to the Qq=12 GeV 2 used in this analysis. The red points indicate the value 
of S p at the average Q 2 of the input data. 



lower edge of the wide range of previous results. Note that this extraction was performed 
using a single model of the deuteron structure and the CD-Bonn N-N potential. No 
estimate of the model dependence due to choice of nuclear models was included. 

Most previous extractions of the neutron structure function treated the data as though 
it were at a fixed Q 2 when in fact there is a strong correlation between x and Q 2 as one 
can see in the left panel of Figure 4. Furthermore, even the simple convolution models 
including only Fermi motion yield a significant Q 2 dependence at large x. This can be 
seen for the calculation of the smearing ratio in right panel of Figure 4 where S p (x) 
is given for different values of Q 2 and, for comparison, S p (x) points at the Q 2 of the 
data are overlaid. Neglecting this x-Q 2 correlation by applying nuclear corrections at 
fixed Q 2 can yield a significant error in the extraction. There are two ways to address 
this correlation. The first is to evolve the model of the deuteron to the Q 2 value of each 
data point, which is essentially the procedure used in the CTEQ6X global extraction 
of parton distributions [7]. The second, which is used in this work is to interpolate all 
of the data to a fixed value of Q 2 . A benefit of this method is that the analysis can be 
performed purely in terms of structure functions, without invoking the parton picture, 
which leads to a particularly direct and accessible analysis scheme. Both of these recent 
approaches yield a similar result, with R np approaching 1/4 (and thus the d(x) /u(x) ratio 
approaching zero) as x becomes very large. 

EXPLORING DEUTERON MODEL-DEPENDENCE 

Given the observation that accounting for <2 2 -dependent contributions is important in 
these analyses [1, 7], we perform a detailed study of the model dependence of the 
extracted neutron structure function when fully and consistently accounting for these 
effects in the comparison. We factor this problem into two parts: 1) the choice of a N-N 
potential and 2) the choice of a nuclear model. We begin by updating the interpolation of 
the global R^ p measurements to fixed Q 2 , this time choosing Q 2 =16 GeV 2 , as this more 



closely matches the kinematics of the large-x measurements, where the Q 2 dependence is 
largest. While this should have no impact on the comparison of the different extractions, 
it should minimize the corrections in the interpolation for all of the results. 

We extract R np , repeating the analysis of Ref. [1] at Qq=\6 GeV 2 . We then take our 
input fit to F2 P along with the extracted F*i n as the starting point to calculate the smearing 
ratios (ratio of the proton or neutron structure function in the deuteron to the free 
structure function) for all of the different models, and use this to determine the change 
in the extracted neutron structure function. Variations in R np can be studied for potential 
and model separately by varying the potential with a single baseline model and varying 
the model with a single baseline potential. These two procedures are depicted in the left 
and right panels of Figure 5 respectively. The model of Ref. [1] is used when comparing 
the different N-N potentials, while the CD-Bonn potential is used for comparison of the 
different calculations of the nuclear effects. We take the full range of the results shown 
in each panel as the one-sigma systematic band for the potential and model dependences 
of the extraction. 
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FIGURE 5. Left panel: R np calculated using various N-N potentials assuming the deuteron model 
described previously [1]. Right panel: R„ p calculated using different models [2, 8, 9, 10, 1] with the 
CD-Bonn N-N potential [5]. 

Note that in this comparison, we use models which include Fermi motion, binding, 
and in some cases off-shell effects and the contributions of nuclear pions. All of these 
represent "ground-up" calculations of the deuteron structure based on input proton and 
neutron structure functions and some N-N potential used to calculate the deuteron mo- 
mentum distribution. Some of the earlier models which yielded exceptionally large re- 
sults for R np at large-.* were based on models which neglected Fermi motion and binding 
and yielded a large effect through significant modification of the nucleon structure in the 
nucleus, e.g. the "scaled EMC effect" result of [3]. While such large nuclear effects are 
possible, they are not implemented in a realistic fashion in this approach. They do not in- 
clude any direct calculation of Fermi motion, are scaled down for the deuteron using an 
extremely large assumed nuclear density density, neglect Q 2 dependence in the nuclear 
effects, and effectively assume S n = S p which is not true when there is a significant dif- 
ference in the x dependence of Fi v and F^. Thus, we do not include models that use an 
explicit "EMC effect" for the deuteron, unless they are included (e.g. via the pion contri- 
butions and off-shell effect) along with Fermi motion and binding fashion [2, 8, 10, 7]. 

We note that the issue of the Q 2 variation of the R^p measurements is important even 
for the models including a large EMC effect in the deuteron. The model used to obtain 
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FIGURE 6. The extracted F^n / F20 with the propagated normalization systematic uncertainty band in red 
and the total deuteron theory uncertainty in blue. The yellow band represents the range of prior extractions 
of this quantity. The arrows on the right once again indicate various values in the x=l limit under certain 
theoretical assumptions as described in Figure 1 . 



the highest set of data points from Figure 1 yields a smaller neutron structure function 
at large x when applied to the Rj p measurements interpolated to a fixed Q 2 value. Thus, 
even if one were to include such models, the impact on the neutron structure function is 
smaller than one would expect based on these earlier comparisons. 

Finally then, the systematic uncertainty bands from the experimental extraction (as 
parameterized in Ref. [1]) and the model and potential dependence of the result are 
combined in quadrature. This is shown in Figure 6 compared to the much larger range 
of prior i^/i^p extractions. Based on these results, it appears that the neutron structure 
function is relatively well known even to large x values, under the assumption that there 
are no "exotic" contributions to the deuteron structure beyond the effects included in 
these models. This conclusion may seem to be at odds with the CTEQ6X analysis [7] 
where, even with the inclusion of large-x data, the uncertainties in the high-jc J-quark 
distributions are extremely large. Figure 7 shows our extracted result for the value 
and uncertainty in the d/u ratio for the proton, where we extract this from F^n/Fzp 
by neglecting contributions from strange and heavier quarks. The left panel shows the 
absolute value of the d/u ratio, and it is clear from this that the J-quark distribution 
is becoming very small. The right panel shows the ratio compared to the reference 
fit, and shows that the relatively small absolute uncertainties on F*i n /Fi v or d/u can 
yield >100% uncertainties on the absolute value of the J-quark distributions. Thus, 
the uncertainty on the J-quark distribution is small relative to the size of the dominant 
w-quark distribution, which means that the data can yield significant constraints when 
comparing to predictions of the d/u ratio at large x. However, the fractional uncertainty 
on the J-quark distribution is large, meaning that observables that are sensitive to the 
J-quark pdf at large x are not well constrained. So the impact of the nuclear models 
is small for some purposes but large in others, and the context is clearly important in 




FIGURE 7. The extracted d/u ratio for the proton with extracted error band. The left panel shows the 
absolute value, while the right panel shows the band relative to the central value. 



determining whether we have sufficiently precise knowledge of the neutron structure 
function. 

There are multiple model-independent extractions planned for the neutron structure 
function at large-x as part of the 12 GeV upgrade plan [11]. These will improve the 
precision of the neutron extraction at high x, but more importantly will be sensitive to 
any physics beyond what is contained in these models, and thus are very sensitive to 
some of the more exotic explanations of the EMC effect. 

CONCLUSIONS 

A recent reanalysis of inclusive DIS data from proton and deuteron scattering, interpo- 
lating to a common value of Q 2 provided a new extraction of F2 n /Fip- The result was 
found to decrease as x becomes large and shows no sign of plateau, sitting near the 
low end of the wide range of previously extracted large-x behaviors. However, some 
previous extractions yielded a neutron structure function that was too large because of 
inconsistent treatment of x-Q 2 correlations. We have extended this analysis to include 
a detailed extraction of the model dependence, using a range of deuteron models and 
potentials. We find that the sensitivity of the extracted F^jFip, and hence the model 
uncertainty, is relatively small, and significantly less than previously believed. While 
these results could be an underestimate if there is a larger than expected "EMC effect" 
in the deuteron, the comparison of this extraction with future model-independent mea- 
surements will be able to measure or significantly limit such modification of the nucleon 
structure in the deuteron. 

ACKNOWLEDGMENTS 

The authors would like to thank Roy Holt for initiating this work and Wally Melnitchouk 
for providing many of the calculations that are crucial in this analysis. This work is 
supported by the U.S. Department of Energy, Office of Nuclear Physics, under contract 



DE-AC02-06CH11357. 



REFERENCES 

1. J. Arrington, F. Coester, R. J. Holt, and T.-S. H. Lee, Journal of Physics G: Nuclear and Particle 
Physics 36, 025005 (2009). 

2. W. Melnitchouk, and A. W. Thomas, Physics Letters B 377, 11 - 17 (1996), ISSN 0370-2693. 

3. L. W. Whitlow, E. M. Riordan, S. Dasu, S. Rock, and A. Bodek, Phys. Lett. B282, 475^182 (1992). 

4. L. L. Frankfurt, and M. I. Strikman, Phys. Rept. 160, 235^127 (1988). 

5. R. Machleidt, Phys. Rev. C 63, 024001 (2001). 

6. J. Gomez, R. G. Arnold, P. E. Bosted, C. C. Chang, A. T. Katramatou, G. G. Petratos, A. A. Rahbar, 
S. E. Rock, A. F. Sill, Z. M. Szalata, A. Bodek, N. Giokaris, D. J. Sherden, B. A. Mecking, and R. M. 
Lombard-Nelsen, Phys. Rev. D 49, 4348-4372 (1994). 

7. A. Accardi, M. E. Christy, C. E. Keppel, W. Melnitchouk, P. Monaghan, J. G. MorfTn, and J. F. 
Owens, Phys. Rev. D 81, 034016 (2010). 

8. S. A. Kulagin, and R. Petti, Nucl. Phys. A765, 126-187 (2006). 

9. A. S. Rinat, M. F. Taragin, and M. Viviani, Nucl. Phys. A784, 25-48 (2007). 

10. Y. Kahn, W. Melnitchouk, and S. A. Kulagin, Phys. Rev. C79, 035205 (2009). 

11. See contributions by G. Petratros, P. Souder, and S. Bueltmann (201 1), this volume. 



